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ABSTRACT 

Aims. We obtained VLT/FLAMES+UVES high-resolution, fibre-fed spectroscopy of five young massive clusters (YMCs) in M83 
(NGC 5236). This forms the basis of a pilot study testing the feasibility of using fibre-fed spectroscopy to measure the velocity dis- 
persions of several clusters simultaneously, in order to determine their dynamical masses. In principle, this reduces the telescope time 
required to obtain a statistically significant sample of dynamical cluster masses. These can be used to assess the long-term surviv- 
ability of YMCs by comparing their dynamical and photometric masses, which are necessary to ascertain the potential evolution of 
YMCs into second-generation globular clusters. 

Methods. We adopted two methods for determining the velocity dispersion of the star clusters: cross-correlating the cluster spec- 
trum with the template spectra and minimising &x 2 value between the cluster spectrum and the broadened template spectra. We also 
considered both red giant and red supergiant template stars. Cluster 805 in M83 (following the notation of Larsen) was chosen as a 
control to test the reliability of the results obtained by this observational method, through a comparison with the results obtained from 
a standard echelle VLT/UVES spectrum obtained by Larsen & Richtler. 

Results. We find no dependence of the velocity dispersions measured for a cluster on the choice of red giant versus red supergiant 
templates, nor on the method adopted. However, we do find that the standard deviation of the results obtained with only one method 
may underestimate the true uncertainty. We measure a velocity dispersion of <X| os = 10.2 ± 1.1 km s for cluster 805 from our fibre- 
fed spectroscopy. This is in excellent agreement with the velocity dispersion of cr [os = 10.6 ± 1.4kms~' determined from the standard 
echelle UVES spectrum of cluster 805. Our FLAMES+UVES velocity dispersion measurement gives A/ V k = (6.6 ± 1.7) x 10 s M Q , 
consistent with previous results. This value of the virial mass is a factor of ~ 3 greater than the cluster's photometric mass, indicating 
a lack of virial equilibrium. However, based on its effective star formation efficiency, the cluster is likely to virialise, and may survive 
for a Hubble time, in the absence of external disruptive forces. Unfortunately, our observations of the other M83 star clusters have 
insufficient signal-to-noise ratios to determine robust cluster velocity dispersions. 

Conclusions. We find that reliable velocity dispersions can be determined from high-resolution, fibre-fed spectroscopy. The advan- 
tages of observing several clusters simultaneously outweighs the difficulty of accurate galaxy background subtraction, providing that 
the targets are chosen to provide sufficient signal-to-noise ratios, and are much brighter than the galaxy background. 
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1. Introduction 

The production of luminous, compact star clusters is character- 
istic of starburst galaxies. However, these clusters have been 
observed in seve ral different ty pes of galaxy, including normal 
spirals (see e.g. iLarsenl 1 19991 and references therein). These 
bright clusters are termed Young Massive Clusters (YMCs), and 
are usually defined as those star clusters that are younger than 
~ 1 Gyr, more massive than ~ 10 s M and compact, with half- 
light radii of a few to a few tens of parsecs. The sizes, lumi- 
nosities and masses of these YMCs are consistent with the prop- 
erties expected of young globular clusters. This led to the sug- 
gestion that YMCs represent globular clusters at an early phase 
of the ir evolution (IWhitmore et alj 119931 : ISchweizer & Seitzeil 
119931) . This would mean that studying YMCs would be tanta- 
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mount to observing globular clusters at the epoch of their birth. 
Since globular clusters are amongst the oldest building blocks in 
our Galaxy, and are benchmarks for stellar and galactic evolu- 
tion, understanding their formation and evolution is fundamen- 
tal in developing our understanding of the fields of large-scale 
star formation, as well as how galaxies build up and evolve. 
Determining whether or not YMCs are proto-globular clusters, 
therefore, is a vital step in gaining insights into these questions. 

The simplest criterion that YMCs need to fulfil in order to 
evolve into globular clusters is that they must be able to sur- 
vive for a large fraction of a Hubble time. Since a cluster re- 
quires sufficient low-mass stars to survive beyond a few Gyr, the 
present-day mass function (PDMF) of a cluster can be used to 
assess its survivability. It may be possible to get a handle on the 
PDMF of a cluster by comparing its dynamical mass to its pho- 
tometric mass, which is the mass predicted by evolutionary syn- 
thesis modelling based on the observed luminosity of the clus- 
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ter and assuming an initial mass function (IMF). If we assume 
that the dynamical mass of a cluster represents its true mass, 
any discrepancies between the two masses can be attributed to 
adopting incorrect assumptions in the photometric mass deter- 
mination, such as the cluster IMF. Even for a very young clus- 
ter, which is not in virial equilibrium, its virial mass can still 
be used to assess the survivability of the cluster. This is done 
by by determining how far out of virial equilibrium the cluster 
is, as quantified by its effe ctive star-formation efficiency (eSFE; 
Goodwin & Bastian 2006). Thus, dynamical mass determina- 
tions can potentially be used to test the scen ario with respect 
to the long-t erm survivability of YMCs, ( e.g. iHo & Filippenkol 
1996a b; see Ide Gri is & Par mentierl 120071 for an overview and 
references therein). 

Since the velocity dispersion of a cluster needs to be mea- 
sured in order to determine the dynamical mass of a clus- 
ter (see Section |3), high-resolution spectroscopy of the clus- 
ter is needed. Large samples of dynamical cluster masses need 
to be obtained to determine whether YMCs might be proto- 
globular clusters, and this can be very expensive in terms of 
telescope time. Observing several clusters at once with a high- 
resolution, fibre-fed spectrograph limits this cost. However, it 
is unclear whether a sufficiently accurate galaxy background 
subtraction can be achieved from only a few fibres placed at a 
different position to the star cluster. In this paper, we explore 
the efficiency of this approach, using FLAMES+UVES to ob- 
serve clusters in M83. Similar attempts at velocity dispersion 
measurements have been made with FLAMES +GIRAFFE (e.g. 
iMieske et al.l l2~008). which uses an order of magnitude more fi- 
bres than FLAMES+UVES. Since large numbers of fibres can 
be dedicated to galaxy background observations, background 
subtraction is more straightforward through the use of a mas- 
ter background spectrum. This, however, is at the expensive 
of at least a factor of two in spectral resolution, and a signifi- 
cantly reduced spectral range. This reduced spectral resolution 
would, however, not permit the measurement of the small clus- 
ter velocity dis persions ant i cipate d for YMCs (~ 7-15 km s ). 
Furthermore, Mies ke et al.l ((2008) observed ultra-compact dwarf 
galaxies (UCDs), rather than clusters. These are much brighter 
relative to the background, and have a much uniform background 
than a spiral like M83. 

M83, also known as NGC 5236, is a prime exampl e of a 
nearby spiral galaxy hosting several YMCs (|L arsen 2004). At a 
distance of 4.5 + 0.3 Mpc (iThim et al.ll2003l) . with its face-on 
inclination and multitude of clusters, M83 is an ideal site to test 
whether fibre-fed spectroscopy can be used to measure the ve- 
locity dispersions of several clusters simultaneously. A further 
advantage of using M83 as the testing ground for this method is 
that the dynamical masses o f two of its cl usters - clusters 502 
and 805 (nomenclature from Larsen 2004) - have already been 
measured from standard high-resolution echelle spectroscopy 
dLarsen & Richtler| [2004). Therefore, these clusters can provide 
a control sample to assess the reliability of the results found from 
fibre-fed spectroscopy. 

M83 is a metal-rich galaxy, wit h a central oxygen a bundance 
of 12 + logO/H = 8.94 ± 0.09 dBresolin et al|2005l). nearly 
twice the solar value of 12 + logO/H = 8.69 dAsplund et al.l 
2004 ). The grand-d esign spiral has a companion, NGC 5253 
dRogstad et al.1 [1974}). with both galaxies hosting intens e star- 
burst activity (ICalzetti et al.1 119991: iTremonti et al.ll200ll) . M83 
has a nuc lear starburst and ongoing star formation within its spi- 
ral arms (Elme green et al.ll998UHarris et al.l200ll) and has a star 
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Fig. 1. Archival FORS 1 image of M83 with the five cluster fibre 
positions and three galaxy background fibre positions indicated. 

formation rat43 of 4 + 1 M yr" 1 . The galax y hosts in excess of 
1000 Wolf-Rayet stars (lHadfield et al.ll2005l) . 

This paper is structured as follows. In Section|2] we present 
our observations of M83 and subsequent data reduction, as well 
as photometry of the clusters and galaxy background regions. 
In Section|3] we determine the velocity dispersion of cluster 805 
and compar e it with the results obtain ed from the standard U VES 
spectrum of lLarsen & Richtferl d2004l) . We discuss the possibility 
of using fibre-fed spectroscopy to obtain velocity dispersions for 
other clusters in Section [4] Finally, we summarise our findings 
in Section[5] 

2. Observations and data reduction 

We obtained fibre-fed spectroscopy of five clusters and three 
regions of galaxy background in M83 with the Fibre Large 
Array Multi Element Spectrograph (FLAMES) linked with the 
Ultraviolet and Visual Echelle Spectrograph (UVES) on the Very 
Large Telescope (VLT) Kueyen Telescope (UT2) in Chile. 

2.1. Imaging of M83 

We retrieved two archival /-band images of M83 taken with the 
Focal Reducer and low dispersion Spectrograph #1 (FORS1) on 
the VLT. These images, taken as part of ESO proposal I.D. 60. A- 
9203(A), were observed on 9th and 11th March 1999 UT for 
180s each. The / filter is centred at 7680A with a 1380A full- 
width at half-maximum (FWHM), similar to the wavelengths 
covered by our FLAMES+UVES spectroscopy. 

The images were obtained to aid in selecting the clusters and 
regions of galaxy background to observe with the eight fibres of 
F LAMES +U V ES . The clusters were chosen from the catalogue 
of lLarsenl d2004l) . First, we selected cluster 805 as a 'control' to 
assess the accuracy of our velocity dis persion measurements, by 
comparing our results with those of lLarsen & Richflerl d2004l) . 

1 This was determined using the equation of Kennicutt] dl998T) and 
adopting the average Ho- flux of the galaxy as 7.1 x lO^'ergs -1 
(Kennic utt, priv. comm.) and a representative extinction of E(B - V) ~ 
0.5 mag dHadfield et al.ll2005h . 
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Table 1. Table of photometry (taken from Larsen 2004), 'background factors' and the average background contribution for the five 
clusters observed. 
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V-I 


814W 


Low 
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contribution (%) 


347 


18.27 


-0.62 


-0.06 


0.40 


19.61 


1.0 


0.7 


0.5 


79 


645 


18.16 


-0.38 


0.23 


0.39 


19.56 


1.2 


0.8 


0.6 


89 


747 


17.36 


-0.65 


0.25 


0.57 


19.80 


1.4 


1.0 


0.7 


76 


805 


16.60 


-0.55 


0.19 


0.65 


17.68 


1.5 


1.1 


0.8 


62 


897 


16.95 


-0.16 


0.51 


0.75 


19.83 


1.7 


1.3 


0.9 


86 



We then short-listed bright candidate clusters that were denoted 
in the catalogue as both iikely/certain cluster' and having been 
successfully f itted with the light-profile fitting program ishape 
dLarsenll 19991) . The final cluster selection was made on the basis 
of successfully positioning the fibres on five clusters and three 
nearby regions of galaxy background with the FLAMES fibre- 
positioning software, prioritising clusters that were flagged as 
isolated with smooth galaxy backgrounds in the catalogue. 

The final select ion comprise d the five clusters 347, 645, 
747, 805 and 897 dLarsenl [2004h and three regions of galaxy 
background, as indicated in Fig. [TJ The galaxy background re- 
gions were observed so that they could be used to subtract the 
galaxy background from the cluster spectra. Therefore, these had 
a range of intensities, similar to the galaxy backgrounds adja- 
cent to the selected clusters. The 'high' fibre was positioned on 
a region with a strong galaxy background, the 'mid' fibre was 
positioned on a mid-intensity galaxy background region, while 
the fibre labelled 'low' was positioned on a region with a faint 
galaxy background. 

The FORS1 images were also used to assess the relative 
fluxes of the galaxy background regions adjacent to the clus- 
ters and within the galaxy background fibres. The images were 
debiased and flat-fielded with standard calibration frames. We 
calibrated the astrometry of the images in Starlink's gai a, using 
the co -ordinates of the clusters given in the catalogue of lLarsenl 
(2004), which were obtained from Hubble Space Telescope 
(HST) imaging. We carried out our aperture photometry on both 
images, also using gaia. We considered the regions of each 
galaxy background fibre as well as the regions of galaxy back- 
ground adjacent to each cluster using several small apertures, 
with radii of both one and two pixels. Taking a mean of these 
galaxy background values, we computed the 'background fac- 
tor' for each cluster and galaxy background fibre pair from the 
relative galaxy background fluxes measured. This background 
factor is the factor by which the galaxy background flux in the 
galaxy background fibre needs to be multiplied in order to (ap- 
proximately) correspond to the galaxy background adjacent to 
the considered cluster. These background factors are contained 
in Table [p a long with photometry of the clusters taken from 
iLarsenl (120041) . 

2.2. Spectroscopy of M83 

We observed the five clusters and three galaxy background re- 
gions of M83 indicated in Fig. QJ over several nights in service 
mode in 2006, using FLAMES +UVES. The average seeing was 
~ 1 arcsec. When linked with UVES, FLAMES is a multi- 
object, high-resolution spectrograph (R = 47 000) with eight 
fibres of 1 arcsec diameter aperture. In order to position the fi- 
bres, each of the eight objects that can be observed at one time 
must have at least a 10.5 arcsec separation. There are two plates 
on the fibre positioner, so that the fibres can be positioned on 



one plate while observations are taking place on the other, thus 
reducing acquisition time. Each plate has a different spectral re- 
sponse and a slightly different wavelength calibration. Although 
UVES is a cross-dispersed echelle spectrograph, with both a blue 
and red arm that can obtain data simultaneously, only the stan- 
dard red UVES set-ups (centred on 5200, 5800 or 8600A) can 
be used with FLAMES. This red arm comprises a mosaic of an 
EEV CCD for the shorter wavelengths and a MIT-LL CCD for 
the longer wavelengths. We observed for a total exposure time 
of 22 048 s, using cross disperser #4, centred at 8600A with 
312 1 mm' 1 , thus obtaining data from ~ 6670 - 10425 A. The 
region ~ 8545 - 8645 A was not observed because it lay within 
the gap between the two CCDs. ThAr arcs were taken for wave- 
length calibration. The resolution of the data, as measured from 
the ThAr arcs, is ~ 9kms~'. The unbinned data have a pixel 
scale of 0.182 arcsec pixeL 1 . 

Since the spectra of M83 were faint, we decided that it was 
more appropriate to reduce the data by hand, rather than to use 
the pipeline-reduced data. The data were reduced using the soft- 
ware package iraf, considering the data from each night, and, 
where relevant, each fibre-positioning plate separately. First, a 
median bias frame for each night was produced and subtracted 
from the data. We then divided each data frame by a master slit 
flat-field for each night. This was produced by taking a median 
of the nine slit flat-fields that had been taken in three different 
slit positions each night, rejecting the six faintest frames, and 
normalising. By rejecting the six faintest frames, the pixels that 
lay in any of the inter-order gaps in the slit flat-fields were re- 
jected, so that most of the chip was illuminated in the master 
slit flat-field. All data lay on the illuminated part of this master 
slit flat-field. The next step was to divide by 'fibre flat-fields', 
which are flat-fields that are obtained through the fibres, and can 
be used to correct for the relative response of each fibre. We di- 
vided by the normalised 2-D all-fibre flat-field. This compared 
favourably with the results achieved by dividing the extracted 
data by the fit to the response of extracted odd- and even-fibre 
flat-fiel ds. All eight fibres of the data were then optimally ex- 
tracted dHorndl 19861) . each with an aperture of 1.456 arcsec (8 
pixels) and wavelength calibrated with the extracted ThAr arcs. 
We determined the heliocentric radia l velocity for each night's 
data using the Starlink program rv d Wallace & Clayton] [ 1997) 
and corrected with these values. The data for each fibre were 
then co-added, and the orders merged. The 1-D data were then 
galaxy-background subtracted. This was done by subtracting an 
average galaxy background spectrum, produced from the two 
galaxy background fibres with sky factors closest to unity (see 
Table [TJ, after multiplication by these sky factors. The relative 
contribution of the galaxy background to the cluster + galaxy 
background spectrum is also included in Table [TJ 

We measured a recessional velocity of v r = 500 + 4 km s , 
after heliocentric radial velocity correction, from the observed 
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central wavelengths of the two Ca n triplet lines visible in our 
data (8498.02 A, 8662.13 A) for cluster 805. These central wave- 
lengths were found by fitting the lines with t he elf (emission-lin e 
fitting) routine in Starlink's dipso package dHowarth et al.l l2004). 
allowing line width, centre and intensity to vary freely. The re- 
cessional velocity implied from the better signal-to-noise ratio 
8498 A line was more heavily weighted, according to the un- 
certainty in the fitted central wavelength estimated by elf. This 
value is in good agreement with other opt ical measurements: 
503 ± llkms" 1 dde Vaucouleurs et al.l|l99ll) . 497 + lSkms -1 , 
445 + 20kms~]_ dFpuque et all Il992l) and 491 + 30kms _1 
dHumason et alj|1956l) . Since no strong lines, including the Ca n 
triplet, could be identified in the spectra of any cluster except 
cluster 805, it was not possible to determine individual cluster re- 
cessional velocities. Therefore, even though small differences in 
radial velocity may exist between the clusters, all cluster spectra 
were velocity corrected by v r = 500 km s . The spectra were 
then normalised, by fitting a polynomial to the continuum across 
the whole spectrum and dividing by this fit. 

2.3. Spectroscopy of standard stars 

With the setup described above (Section [2. 21 . we also observed 
a telluric star, CD-32°9927, and three Galactic red giants - 
HD 159821 (Kl III), HD 123833 (K2 III) and HD 159881 
(K5 III). We chose to observe giants rather than supergiants, be- 
cause these four calibration stars could be observed in only two 
fields, reducing the overheads of the observations. For each of 
these standard stars, one sky fibre was placed as close to the ob- 
ject as possible. 

We chose to use the pipeline-reduced data for these bright 
standard stars, since these included a better extraction of the 
8662 A Ca n triplet lines than the manually-reduced data. This 
line lies very close to the edge of the MIT-LL CCD for the tem- 
plate stars. Due to their higher recessional velocity, this is not 
a problem for the M83 clusters. We sky subtracted the spec- 
tra of these bright standards by simply subtracting the relevant 
sky fibre spectrum. The data were heliocentric radial velocity- 
corrected using the values computed by rv. Each red giant was 
then corrected for its recessional velocity using the observed 
central wavelengths of the two visible Ca n triplet lines, as mea- 
sured with elf. The spectra were then normalised. 

An example spectrum of a velocity-corrected, normalised 
red giant template star, HD 159821 (Kl III) is shown in Fig. [2] 
for a reduced wavelength coverage. The principal telluric fea- 
tures in this spectral range and the two visible Ca n triplet lines 
are marked. An example cluster spectrum over suitable spectral 
regions is shown in Fig. [3] fitted with a red giant template star. 

3. Virial mass 

The virial mass of a cluster, M vn - is given by the equation: 
M- 

vir G ' 

where <xi os is the line-of-sight velocity dispersion of the cluster 
and ru is the projected half-light radius of the cluster. This equa- 
tion follows from the virial theorem, with rj = 9.75 for a wide 
range of light profiles, assuming that the cluster is virialised, 
gravitationally bound and spherically symmetrical, that the ve- 
locity dispersion of the cluster is isotropic and that all of t he stars 
contained in the cluster are single stars of equal mass dSpitzerl 



119871) . However, recent studies have shown that 77 is not a con- 
stant, and can vary with time, depending on, for example, the 
degree of mass segregation and the binary fraction of the cluster 
(Boil v et all [20051 iFleck et al.1 120061: iKouwenhoven & de Grijsl 
2008). So that our results are comparable with other dynamical 
mass determination of YMCs, we adopt 77 = 9.75, as has been 
done by other authors. 

3.1. Velocity dispersion 

Clearly, it is necessary to measure the line-of-sight velocity dis- 
persion of the cluster in order to determine the cluster's virial 
mass. This is achieved by comparing the widths of the Doppler- 
broadened lines of the cluster with a template star. The cluster 
spectrum is compared to either red giant or red supergiant tem- 
plate star spectra, because the lines of these stars are broadened 
by only a few km s" 1 by macro-turbulence in their atmospheres 
dGrav & Tonerll986h . It is not appropriate to use earlier type su- 
pergiants and main sequence stars, as effects such as rotational 
broadening, macro-turbulence and micro-turbulence broaden the 
lines of these stars by amounts substantially greater than the an- 
ticipated cluster velocity dispersions. Therefore, only spectral 
regions redwards of ~ 5000 A, which are dominated by light aris- 
ing from cool su pergiants of spectral typ e F-M, should be con- 
sidered (see e.g. lHo & F ilippenko Jl996al) . Red supergiant tem- 
plates are preferable, as young clusters do not contain red giants. 
However, red giants can be used, although these have lower val- 
ues of macro-turbulent broadening than more luminous giants. 
This could mean that comparing the spectrum of a cluster that is 
dominated by red supergiants with a red giant template spectrum 
may syste matically overestim ate the cluster velocity dispersion. 
However, lLarsen et al.l (|2004) found no systematic differences, 
within the uncertainties, in the velocity dispersions determined 
for a cluster when cross-correlating with templates of luminosity 
classes I, II and III. 

Here, we adopted two methods to determine the line-of-sight 
velocity dispersion. The first involved minimising ay 2 value be- 
tween the cluster spectrum and broadened template spectra. The 
second involved cross-correlating the cluster spectrum with each 
template spectrum. These two methods are discussed in more de- 
tail in Sections IXTTTl and |3~L2l 

To make a direct comparison between the results that can 
be obtained with fibre-fed spectroscopy and those from a stan- 
dard echelle spectrum, we also determined velocity dispersions 
using both methods with the non-fibre-fed UVES spectrum of 
cluster 805 from lLarsen & Richtlerl d2004l) . We degraded the 
UVES spectrum by re-binning the data to the same linear dis- 
persion scale as our own. To consider the effect of red super- 
giant over red giant templates we also included two red su- 
pergiant templ ate stars from the U VES Paranal Observatory 
Project (POP; iBagnulo et ai1l2003l) : HD 206778 (K2 lb) and 
HD 12642 (K5 lab). These were also re-binned to the same lower 
linear dispersion scale as our own data. 

3.1 .1 . x 1 minimisation 

We measured the line-of-sight velocity dispersion of cluster 805 
by minimising a y 1 value between the cluster spectrum and 
each broadened template spectrum in turn. All spectra were nor- 
malised and then the continuum was subtracted. The broaden- 
ing was achieved by convolving a template spectrum with a 
Gaussian of <x equal to the desired velocity broadening. To find 
the best fit, the ^-minimisation program considered a range of 
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Fig. 2. Pipeline-reduced spectrum of the Kl giant HD 159821, normalised and corrected to rest wavelength, over the wavelength 
range 6750 - 9150 A. The two lines of the Ca n triplet that are contained in our data (8498.02 A, 8662.13 A) are marked. The other 
line of the Ca n triplet (8542.09 A) lies on the gap between the two CCDs, which is also indicated. The major telluric features in 
this spectral range are also marked. 



broadening values, over selected spectral regions. The criteria 
for selection were that the spectral regions showed visible simi- 
larity in absorption lines between the cluster and template spec- 
tra and that they did not contain any telluric lines, which could 
provide an artificial match between the spectra. These spectral 
regions were: 7335 - 7370A, 7403 - 7445 A, 7460 - 7513A, 
8348-8385A, 8465-8502Aand 8659-8712A. Rather than con- 
sidering these small regions of spectra individually, we ran the 
model over the whole spectrum, but only considered the chosen 
spectral regions when computing the y 2 value. Non-linear steps 
of broadening were considered in the ^-minimisation program, 
with closer steps as the minimum was approached. An exam- 
ple fit of a broadened template star spectrum to the cluster 805 
spectrum over these spectral regions is shown in Fig. [3] 

Since the Ca n triplet lines may be saturated in the template 
spectra, including the m in the fit may overestimate the cluster ve - 
locity dis persion (e.g.lWalcher et al.ll2005t iMartini & Hol l2004). 
However, iMengel et al.l ( 20021) found no disparity between the 
velocity dispersion results computed by x 2 minimisation for the 
strongest component of the Ca n triplet and other individual ab- 
sorption features for clusters in NGC 4038/4039. We decided to 
consider the potential impact of this effect by excluding these 
lines from the fit. Therefore, we considered three cases, with 
each case including slightly different spectral regions in the fit. 
The first case used the spectral regions described above, includ- 
ing both visible Ca n triplet lines. Thus, it was called 'Both CaT' . 
We also considered the case where the 8498A Ca n line was 



masked out of the fit, still including the 8662A Ca n line, since 
the 8498A was not in the UVES spectrum of lLarsen & Richtierl 
d2004t) . This case was called '8662 A' . In the final case, both Ca n 
lines were masked (called 'No CaT'). The results for all of these 
cases are given in Table [2] 
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Fig. 3. Best fit of broadening for the K2 III template HD 123833 (blue line) to our normalised FLAMES+UVES cluster 805 
spectrum (black line) for the 'Both CaT' case of 10.9 km s . The spectral regions included in the fit are shown where th e red line is 
at 1 .5. Also included in this figure is the standard echelle UVES spectrum of cluster 805 from lLarsen & Richtlerl (120041) (grey line), 
offset at -1. See the online journal for a colour image. 



The results obtained from this method may be sensitive to 
low signal- t o-nois e ratios. Taking Ax 2 = 1 as the Icr error 
dPress et alj Il992l) . we estimate that the uncertainties on the 
FLAMES+UVES velocity dispersions are ~ 2kms~ 1 , while the 
uncertainties on the UVES velocity dispersions are ~ 1 kms" . 



3.1.2. Cross-correlation method 

In this second method, we determine the clust er velocity dis- 
persio n using the cross-correlation technique of iTonrv & Davisl 
d 19791) . All spectra were normalised and then the continuum was 
subtracted. We cross-correlated the cluster spectrum with each 
template spectrum in turn, over the spectral regions described in 
Section [3. 1.1 1 as 'No CaT', excluding the Ca n triplet lines. We 
did not consider the other two cases, since broadly consistent 
results for all three cases were obtained for the x 2 minimisation. 

The FWHM of the resulting cross correlation function (CCF) 
relates to the velocity dispersion of the cluster. This relationship 
was calibrated by considering the template star used. The tem- 
plate star was broadened and cross-correlated with the original, 
unbroadened template. By considering several values of broad- 
ening, and measuring the FWHM of each CCF, the almost linear 
relationship between the FWHM of the CCF and broadening was 
empirically calibrated to an absolute scale, from which the ve- 



locity dispersion of the cluster was read. This was repeated for 
each template, because every star has a different calibration. 

This method is less sensitive to spectral-type matching and 
low signal-to-noise than the x 2 minimisation. However, it suf- 
fers from complications associated with the subjectivity of fit- 
ting weak, non-Gaussian CCFs. Based on a 5 % error in fitting 
the FWHMs of the CCFs, each measurement has an uncertainty 
of ~ lkms . Note that it is not the intrinsic FWHM that is 
important, but the relative FWHMs. 



3.2. Comparison of results and techniques 

Firstly, it is apparent from Table [2] that there is no systematic 
difference in the velocity dispersion measured when comparing 
a cluster spectrum with red giant or red supergiant templates. 
This is true for both the^ 2 minimisation technique and the cross- 
correlation. For this reason, we will consider the mean of all five 
templates when discussing results below. 

Secondly, we do not observe a systematic increase in the ve- 
locity dispersion measured when the Ca n lines were included 
in the^ 2 minimisation. This may be due to the large wavelength 
range considered in the fit, thus reducing the weight applied to 
the strong Ca n triplet lines. This would also explain the poor fit- 
ting of these lines. Due to this similarity in results, we consider 
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Table 2. Velocity dispersions in kms" 1 determined by comparing our FLAMES +UVES cluster 805 spectrum and the standard 
UVES cluster 805 spectrum of lLarsen & Richtierl d2004l) with five red giant (RG) and red supergiant (RSG) template spectra. 





X 1 minimisation 


Cross-correlation 


Template 


UVES 




FLAMES +UVES 


UVES 


FLAMES+UVES 




'No CaT'* 


'8662A'* 


'No CaT'* 


'8662A*' 


'Both CaT'* 


'No CaT'* 


'No CaT'* 


Kl III 


10.3 


10.4 


9.5 


9.8 


9.9 


9.8 


9.1 


K2 III 


10.8 


10.9 


10.2 


10.6 


10.9 


9.8 


10.0 


K5 III 


11.8 


11.9 


10.8 


11.1 


10.9 


8.4 


8.4 


K2Ib 


13.2 


13.2 


12.5 


12.7 


12.6 


9.8 


10.4 


K5 lab 


11.7 


11.7 


10.8 


11.1 


11.0 


10.0 


10.2 


Mean of RGs 1 


10.9 (0.8) 


11.0 (0.8) 


10.1 (0.7) 


10.5 (0.7) 


10.6 (0.6) 


9.3 (0.8) 


9.2 (0.8) 


MeanofRSGs 1 


12.4 (1.1) 


12.5 (1.0) 


11.7 (1.2) 


11.9 (1.1) 


11.8 (1.1) 


9.9 (0.1) 


10.3 (0.1) 


Mean of all 1 


11.5 (1.1) 


11.6 (1.1) 


10.8 (1.1) 


11.1 (1.1) 


11.1 (0.9) 


9.6 (0.7) 


9.6 (0.8) 



*This indicates the spectral region considered. 'Both CaT' includes all selected spectral regions (see text), '8662A' excludes the 8498A Ca n 

line and 'No CaT' excludes both visible Ca n lines. 

1 The standard deviation on each mean is presented in parentheses 



the results obtained for all five template stars with the 'No CaT' 
case of spectral regions for both techniques. The mean and stan- 
dard deviation of these results give <xi os = 10.2 + 1.1 km s -1 . 

We find that, within the uncertainties the velocity disper- 
sions found from the FLAMES+UVES data are consistent with 
those measured from the standard UVES data, for both meth- 
ods and all spectral regions considered. If we also take a mean 
of the ten measurements obtained for the two techniques with 
the 'No CaT' case of spectral regions for the standard UVES 
data, we find <x los = 10.6 + 1.4 km s , again with the uncer- 
tainty representing the standard deviation of the results. This is 
in excellent agreement with our FLAMES+UVES measurement 
of cr los = 10.2 + 1.1 kms -1 . Therefore, we conclude that our 
results agree with those obtained from standard echelle spec- 
troscopy, and that fibre-fed spectroscopy is practicable for the 
purpose of determining cluster velocity dispersions. 

Our FLAMES+UVES velocity di spersion result is also in 
agreement with the result measur edbyE arsen &RichtleH (12001 
of cr los = 8.1 + 0.2km s , at the 2<x level. This was determined 
by cross-correlating their UVES spectrum with red supergiant 
templates. They considered several spectral regions separately: 
4520 - 4850A, 73 10 - 7570A, 77 10 - 8070A and 8680 - 8880A, 
as well as considering two cases including multiple regions: the 
region 6770 - 6850A with 6940 - 71 10A and the region 6770 - 
6850A with both 6940 - 71 10A and 7310 - 7570A. 

We only obtained results for cluster 805. The other clusters, 
which are fainter and have higher relative background contribu- 
tions, have lower signal-to-noise ratios. We measure a contin- 
uum signal-to-noise ratio per resolution element (« 0.1 A) for 
cluster 805 of ~ 7 at 7000 A and ~ 2 at 9000 A compared with 
a signal-to-noise ratio per resolution element for cluster 897, the 
next brightest cluster in the /-band, of ~ 6 at 7000 A and < 1 at 
9000 A. For these other clusters, the noise is such that no spec- 
tral regions could be found for which the absorption lines present 
in the template star could be identified in the cluster spectrum. 
Computing the x 2 minimisation over the regions adopted for 
cluster 805 led to a broadening that was an order of magnitude 
too large. Cross-correlation over these regions did not produce 
a peak in the CCF In the case of cross-correlation, the lack of 
results cannot be due to differences in radial velocities between 
these faint clusters. Indeed, had a CCF peak been produced, the 
radial velocities of these clusters could have been measured from 
its position. 



The rapid decline in the signal-to-noise ratio in the spectra 
of all clusters, including cluster 805, redwards of the Ca n meant 
that we did not include these regions in the velocity dispersion 
determinations. If the galaxy background was not subtracted 
from the spectrum of cluster 805, and the sky emission lines 
were masked out of the fit, the velocity dispersion of the cluster 
could not be measured. This is probably due to the fact that the 
lines were significantly diluted by the background. Therefore, 
this could also not be used to measure the velocity dispersion for 
any of the other four clusters. 

Clearly, in order to measure reliable star cluster velocity dis- 
persions, good signal-to-noise ratios in the cluster spectra are 
paramount. We obtained consistent results from cluster 805 with 
a signal-to-noise ratio of > 3 in the spectral regions considered, 
but no results from the data of the other clusters, with lower 
signal-to-noise ratios. Therefore, a signal-to-noise ratio of ~ 3 
seems likely to represent a lower limit to the signal-to-noise ra- 
tio required to obtain reliable results. To investigate this, we pro- 
duced a simple model cluster by adding the library K2 lb and 
K5 lab spectra, re-normalising the spectrum and broadening it 
by 10 kms -1 . A y 2 minimisation was carried out between this 
model spectrum and all five template stars for the 'Both CaT' 
case of spectral region, in an attempt to retrieve this value of 
the velocity dispersion. The signal-to-noise ratio of the model 
cluster spectrum was degraded by artificially adding Gaussian 
noise and the^- 2 minimisation was recomputed for the new nois- 
ier model spectrum. This was repeated for a range of model clus- 
ter spectrum signal-to-noise ratios. 

The results are shown in Fig. |4] The mean of the five values 
measured for each signal-to-noise ratio are also plotted, and an 
uncertainty on this mean equivalent to the standard deviation is 
represented by the dashed lines. The measured velocity disper- 
sion increases with decreasing signal-to- noise ratio of the clus - 
ter, in agreement with the results found by Marti ni & Hoi (12004). 
The plot also shows that, in this case, a cluster signal-to-noise 
of > 2.1 is within 20% of the true broadening, and the lower 
limit on the signal-to-noise ratio for which the standard devia- 
tion of the results describes the true error is just over 3. The pre- 
cise signal-to-noise ratio at which these minima occur will de- 
pend of the cluster and on the nature of the template mismatch, 
which should be less stark than that illustrated here. This plot 
indicates that the velocity dispersion measured for cluster 805 
is robust, although, as already stated, the standard deviation of 
the results from the x 1 minimisation may underestimate the true 
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Fig. 4. Plot of velocity dispersion measured from ax 2 minimisa- 
tion between the five template stars and a simple model cluster 
spectrum with <x los = lOkms -1 in the 'No CaT' case against 
the signal-to-noise ratio (SNR) of the model cluster. The mean 
velocity dispersion measured is also plotted. The dashed lines 
represent the uncertainty, as measured by the standard deviation 
if the values, on this mean. The dotted line shows the true veloc- 
ity dispersion of the model cluster. See the online journal for a 
colour image. 

Table 3. Parameters for cluster 805 taken from the literature and 
determined in this work. 



Parameter 


Value 


Reference 


Distance (Mpc) 


4.5 ±0.3 


Thim et al. (2003) 


My (mag) 


-12.2 ±0.4 


Larsen & Richtler (2004) 


Age (Myr) 


13±6 


Larsen & Richtler (2004) 


M phot (xl0 5 M o )* 


1.9 ±1.4 


Larsen & Richtler (2004) 


rm (pc) 


2.8 ±0.4 


Larsen & Richtler (2004) 


o"i os (kms -1 ) 


10.2 ±1.1 


This work 


M vlr (x 10 5 M Q ) 


6.6 ± 1.7 


This work 


L v /M yk (Lq/Mq) 


10±5 


This work 



* assuming Kroupa (2002) IMF. 



uncertainty. The larger uncertainty measured from the standard 
deviation of all five templates with both methods should remain 
a reliable estimate of uncertainty. 

3.3. Implications for cluster 805 

Adopting the cluster parameters determined for cluster 805 con- 
tained in Table [3] and using the velocity dispersion measured 
here, we find M vir = (6.6 + 1.7) x 10 5 M s and a V-band 
luminosity-to-dynamical mass ratio Ly/Mvir = (1O±5)L /M O . 
These results a gree, within the u ncertainties, with the values 
determined by Larsen & Richtlerl d2004l) of M vil = (4.2 + 
0.67) x 10 5 M Q and L v /M vir = (15 + 6) L o /M . Both of 
these dynamical mass measu rements are la rger than the pho- 
tometric mass predicted for a iKroupal (120021) IMF of M p h t = 



(1.9 ± 1.4) x 10 5 M B for the cluster age of 13 ± 6 Myr 
(lLarsen & Richtlerl |2004|) . indicating that cluster 805 may not 
be virialised. However, the photometric and dynamical masses 
agree at just over the 2<x level. 

Our value of light-to-dynamical-mass ratio suggests that the 
eSFE of cluster 805 lies between ~ 30 % and ~ 50 %, favouring 
the value of ~ 40 %. lGoodwin & Bastianl d2006l) estimate that a 
cluster with an eSFE ~ 30% will lose ~ 80% of its mass within 
~ 40 Myr and disperse. However, a cluster with an eSFE > 40% 
will virialise, although it will lose up to ~ 60% of it mass, and 
could survive for a Hubble time in the absence of external dis- 
ruptive forces. 



4. Discussion: using high-resolution, fibre-fed 
spectroscopy to determine velocity dispersions 

Using high-resolution, fibre-fed spectroscopy to determine ve- 
locity dispersions would have a major impact on producing a 
large sample of dynamical cluster mass determinations, by mak- 
ing more efficient use of telescope time. By having a large sam- 
ple of dynamical cluster masses, it should become possible to 
assess how likely YMCs are to survive for a large fraction of a 
Hubble time and evolve into globular cluster-type objects. 

The primary concern over this method is the question of how 
successfully the galaxy background could be subtracted from 
fibre spectroscopy of faint clusters. Of course, sky subtraction 
is not a problem with bright stars, despite a small amount of 
flux 'bleeding' from the bright star fibres into the sky fibres. For 
the fainter objects, galaxy background subtraction is more chal- 
lenging. Simply scaling and subtracting the galaxy background 
fibre with the most similar intensity background to the clus- 
ter being considered introduces more noise into the spectrum. 
Taking an average of two or three scaled galaxy background fi- 
bres minimises this, although still introduces noise. Subtracting 
an average of two galaxy background fibres decreases the signal- 
to-noise ratio of the cluster + background spectrum by a factor 
of ~ 2.5 for cluster 805, the brightest cluster, and a factor of 
~ 4 for both cluster 347, the faintest cluster, and cluster 645, 
which has the lowest contrast with the galaxy background. It 
seems that scaling the galaxy background in this way is a rea- 
sonable approximation to the local galaxy background. All three 
galaxy background fibres are very similar, considering the level 
of noise, in spite of spatial variations in galaxy background, 
due, for example, to differing stellar populations and extinc- 
tion. Nevertheless, some sky lines are imperfectly subtracted in 
the cluster 805 spectrum, although these are within the level of 
the noise. A reduced^- 2 between the low- and mid-background 
spectra and the scaled, high-background spectrum yields values 
of x 1 = 0.9 and x 1 = 0.7, respectively. While adopting the 
noise in the spectra as the uncertainty in the fit seems to overes- 
timate the true uncertainty, the relative similarity in the galaxy 
background spectra, given the large degree of noise, is demon- 
strated. An alternative method would be to fit the continuum of 
the scaled, extracted galaxy background fibre and subtract this fit 
from the cluster spectrum, manually removing the sky lines from 
the cluster spectrum. However, the extra complications involved 
in this process were not warranted for our data, since an appro- 
priate galaxy background was achieved from averaged galaxy 
background fibres. 

Also, the s ignal-to-noise rat i o was much higher in the UVES 
spectrum of lLarsen & Richtlerl (120041) . for which a shorter in- 
tegration time of 7500 s was used. This would appear to indi- 
cate that long-slit echelle spectroscopy is a more efficient use 
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of telescope time, especially since two clusters could poten- 
tially be observed in the same slit. However, finding cases where 
two clusters can be observed simultaneously can be difficult, 
due to the relatively short slit length of many echelle sepctro- 
graphs (e.g. 1 1 arcsec for UVES). The signal-to-noise ratio in 
the galaxy-subtracted UVES+FLAMES cluster spectra was sub- 
stantially lower than anticipated from the ESO exposure time 
calculator. This predicted a signal-to-noise ratio in excess of 
30 at ~ 8000 A, neglecting the high galaxy background. Even 
considering this additional noise does not account for the much 
lower observed cluster signal-to-noise ratios. The most likely ex- 
planation for the greatly reduced signal-to-noise ratio is is that 
the fibres were positioned off-centre. However, this needs to be 
investigated further. 

As discussed in Section l3~2l the signal-to-noise ratio of clus- 
ter 805 is a good representation of the lowest signal-to-noise ra- 
tio required to measure robust velocity dispersions. Therefore, 
making a conservative assumption that 50% of the flux of clus- 
ter 805 was missed in our spectrum indicates that results could 
be obtained for a cluster half as bright as cluster 805 (i.e., 
mi « 16.7 mag). Approximately 18 of the 159 YMCs i n M83, 
and 10 9 of all 1358 YMCs, contained in the catalogue of lLarsen] 
(2004) fulfil this criterion. 

Another problem with this method, specific to 
FLAMES+UVES, is the limited choice of spectral region, 
as dictated by the need to use a standard UVES setup. This 
meant that one of the strong lines of the Ca n triplet was lost. 
Furthermore, much of the data in the near-infrared cannot be 
used to determine velocity dispersions, due to telluric features 
and the lower CCD response in this spectral region. Clusters 
also need to be selected carefully, with M v i r > 10 5 M o in order 
to have high enough velocity dispersions to be measured. For 
a 10-Myr-old cluster, this corres ponds to My = -11.5 mag 
(starburst99: iLeitherer et aD,ll999l) . or my ~ 16.8 mag for a 
distance of 4 Mpc and E(B -V) ~ 0.1 mag. 

While the drawbacks of this method must still be borne in 
mind, we have shown here that reliable velocity dispersions can 
be measured from fibre-fed spectroscopy. Therefore, this method 
can potentially be used to obtain a large sample of dynamical 
masses, which can be used to address the issue of the nature of 
YMCs as proto-globular clusters. 



5. Summary 

We determined a velocity dispersion of cluster 805 in M83 of 
cr los = 10.2 + 1.1 km s" 1 from VLT/FL AMES +UVES spec- 
troscopy. This agrees well with the results obtained here for 
the standard echel l e VLT /UVES spectrum of cluster 805 from 
lLarsen & Richtlerl d2004h of cr los = 10.6 + 1.4kms^. This 
shows that high-resolution, fibre-fed spectroscopy is a practi- 
cable method with which to measure the velocity dispersions 
of several clusters simultaneously. This allows a more effi- 
cient use of telescope observing time, providing that sufficiently 
high signal-to-noise ratios are achieved. Therefore, high-number 
statistics of YMC dynamical masses, vital for assessing the role 
of YMCs as proto-globular clusters can be more easily obtained. 

We find no systematic differences in the velocity dispersions 
measured depending on whether red giant or red supergiant tem- 
plates were considered. Nor do we find any differences in the 
velocity dispersions determined from the x 1 technique when the 
Ca n triplet lines are included in the fit. While the results pro- 
duced from both a x 1 minimisation and from cross-correlation 
are consistent, there are indications that adopting the standard 



deviation of the results determined from only one technique may 
underestimate the true uncertainties. 

Adopting our value of velocity dispersion and all other 
cluster parameters from lLarsen & Richtlerl (|2004|) . we measure 
a virial mass of M vir = (6.6 + 1.7) x 10 5 M Q for clus- 
ter 805, somewhat larger than its photometric mass of M p h t = 
(1.9 + 1.4) x 1O 5 M , for a iKroupal d2002l) IMF and an age 
of 13 ± 6 Myr. Co mbining this with the eSFE predictions of 
iGood win & Bastianj d2006l) . we conclude that the cluster will 
probably virialise and could survive for a Hubble time in the 
absence of external disruptive forces. 
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